Scattering of plasmons at the intersection of two nanotubes: Implications for 

tunnelling 
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We study theoretically the plasmon scattering at the intersection of two metallic carbon nan- 
otubes. We demonstrate that for a small angle of crossing, 8 < 1, the transmission coefficient is 
an oscillatory function of X/9, where X is the interaction parameter of the Luttinger liquid in an 
individual nanotube. We calculate the tunnel density of states, u(lj,x), as a function of energy, u>, 
and distance, x, from the intersection. In contrast to a single nanotube, we find that, in the geome- 
try of crossed nanotubes, conventional "rapid" oscillations in u(u>,x) due to the plasmon scattering 
acquire an aperiodic "slow-breathing" envelope which has X/9 nodes. 

PACS numbers: 71.10.Pm, 73.40.Gk,72.15.Nj 
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Introduction. By now, observation of Luttinger liquid 
in ID systems has been reported for single- walled car- 
bon nanotubes [l j. 2, & 0, d, @] and GaAs-based semi- 
conductor wires 7j. Conclusions about Luttinger liq- 
uid behavior have been drawn from analysis of the data, 
which can be divided into two groups: (i) power-law, 
oc (max{V, T|)" behavior of tunnel or source-drain con- 
ductance [l|, 2, E3, EL H> @| ) where parameter a is the mea- 
sure of deviation from the Fermi liquid behavior, and 
(ii) momentum-resolved tunnelling in a parallel magnetic 
field 0. 

On the conceptual level, the difference between the 
techniques (i) and (ii) is that (i) probes a single- 
point Green function, Q(x,x,u>), while (ii), by mapping 
J dx J dx'Q(x, x', to) exp[— iq B (x — a/)], with q B propor- 
tional to applied field, yields information about two-point 
Green function, and thus is more informative. 

With regard to quantitative determination of the Lut- 
tinger liquid parameter, g, which is related to a as Q 
a = (g^ 1 + g — 2)/8, it is desirable to identify an effect, 
which would depend on g stronger than a power law. 
An example of such an effect was given by Ussishkin and 
Glazman in Ref. 0, where, due to electron backscattering, 
g appears in the argument of sine; this sine describes the 
amplitude modulation of the probe-induced Friedel oscil- 
lations [lfj, oc cos(2fc F cc) in the local density of states; k F 
is the Fermi momentum. 

In the present paper we demonstrate that the geom- 
etry of the crossed nanotubes (see Fig. 1) offers a qual- 
itatively new manifestation of the Luttinger liquid be- 
havior. In particular, the oscillatory dependence on g, 
similar to that in Ref. 0, emerges in the geometry of 
crossed nanotubes even without electron backscattering 
0, [ll|, Qjl, [l]|. More precisely, we show that, in this 
geometry, the envelope, "breathing" with g, modulates 
not cos(2fc F a;) oscillations, but much slower oscillations 
resulting from the plasmon backscattering. 

There is an important difference between scattering 
of plasmons and electrons: for an obstacle bigger than 
hp 1 electron scattering is exponentially suppressed, while 




FIG. I: (Color online) Intersecting nanotubes, separated by 
a distance, d; the angle of intersection is 6. Directions of 
incident, reflected, r, transmitted, t, and deflected, dz,, da, 
plasmon waves are shown with solid red arrows. Dashed ar- 
rows illustrate two contributions to the reflected wave. 



plasmon scattering is efficient as long as the size of the 
obstacle does not exceed the plasmon wavelength. This 
scattering gives rise to the oscillations of the local density 
of states 5v(u),x) oc cos(2ujx/v f ) where v F is the Fermi 
velocity. It is these oscillations that acquire a breathing 
envelop in the geometry of crossed nanotubes, Fig. 2. 
Our main finding is that, with regard to this modula- 
tion, making the crossing angle 9 small, effectively en- 
hances the Luttinger liquid parameter. To describe this 
enhancement quantitatively, we first consider an auxil- 
iary problem of plasmon scattering at the intersection 
and later utilize it for the calculation of Sv(uj,x). 
Plasmon scattering at the intersection. Assume that d 
is the minimal distance between the nanotubes. Even in 
the absence of electron tunnelling, a plasmon propagat- 
ing towards x — in the nanotube I can: (i) pass x = 
(transmission); (ii) excite a plasmon in the nanotube 2, 
which propagates away from the intersection x — either 
to the left or to the right (deflection); (Hi) get reflected. 
Incorporating the plasmon scattering into the formalism 
of the Luttinger liquid gives rise to the breathing enve- 
lope in Fig. 2. The underlying reason is that the interac- 
tion between the tubes strengthens towards intersection 
flij ]. This leads to the x- dependent splitting of velocities 
in each tube. The resulting x-dependent phase accumu- 
lation near the intersection translates into nontrivial de- 
pendence of 5v(u>,x). Moreover, the phase accumulation 
increases rapidly with decreasing angle 9, thus simulating 
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FIG. 2: (Color online) Oscillating corrections to the tun- 
nelling conductance is plotted from Eq. (|15[) versus dimen- 
sionless bias w = Vx/s, where x is the distance from the 
intersection. Periodic oscillations (black) are modulated by 
"breathing" envelope (red), with "period" determined by X/6, 
where A is the interaction parameter; (d) illustrates suppres- 
sion of oscillations at finite separation, d — 0.01a;, between 
the nanotubes. 



the enhancement of the Luttinger parameter. 
Collective modes of intersecting nanotubes. As a result of 
long-range interaction, e 2 / 9 1^ ) , where u{x) 

is the displacement of the electron position from the equi- 
librium, the plasmon spectrum of an individual tube is 
^(l) — Qs(q) W1 th velocity Q s — v F (l + Aln(qr)) 1 / 2 . 
Here r is the nanotube radius and A = 8e 2 /(Trhv F ) is 
the interaction constant 0. Following Ref. [H, we ne- 
glect the relative change of ln(qr). At a given frequency, 
lu, displacement u(x) is the eigenmode, D{u(x)} = 
(uj 2 /v 2 ) u(x), of the operator 
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For two crossed nanotubes, the eigenmodes are described 
by the system of two coupled equations, 



{lo 2 /v 2 f ) u h2 (x) = D{u h2 (x)} + F{t*2,i(a;)}, (2) 



Hf} 



dx 



df(x) 



dyJ-oo^/d 2 + x 2 + y 2 - 2xycos9 dx 



(3) 



The operator F has a meaning of longitudinal force cre- 
ated by the density fluctuation, du\{x) / dx, in the nan- 
otube 1, at point y of the nanotube 2. The scattering 
problem corresponds to the solution of Eqs. @ which 
has the following asymptotes at large x and y: 
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Born approximation. For small A, the elements of scat- 
tering matrix can be found in the Born approximation in 
momentum space. To the first order in A, only da and 
c?l are non-zero. They are given by matrix elements of 



the operator F, Eq. ([3]), namely, dn(k) = (i/2k)Fk^ and 
dh(k) = —(i/2k)F-k,k- Analytical expression for F p , q is 



F p , q = 2nX 



pqe 



(p 2 +q 2 -2pqcos9) 1/2 



(p 2 + q 2 — 2pq cos 
This leads to the final result for deflection coefficients 



(5) 



dn = iirX- 



2 sin | 



-iirX- 



2 cos I 
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An apparent consequence of Eq. (|6|) is that deflection is 
exponentially small when the plasmon wavelength is <C d. 
Less obvious is that, for kd < 1 and small 0, coefficients 
da and d^ can differ exponentially. This is because the 
exponent, exp(— 2kd/&), in dh can be small if kd is small. 
Noteworthy, in the long-wavelength limit, kd <C 0, we 
still have d^/dh ~ 1/6 » 1. The underlying reason is 
that dh corresponds to the wave which travels almost 
in the opposite direction to the incident wave, while da 
travels almost along the incident wave. From Eq. ([6|) we 
conclude that the Born approximation applies at A <C 9. 

The reflection coefficient, r(k), in the second Born ap- 
proximation, is expressed via the matrix elements Eq. ([5]) 
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The integral in Eq. is the sum, (27rA) 2 fc(/ 1 + il 2 ), of 
contributions from the poles p = ±fc and the principal 
value, which can be cast in the form 



h=P 



dpp 2 e ~inre(P +k 2 -2pkco S 9) 
p 2 - k 2 ^ p 2 _|_ £2)2 _ 4 p 2j e 2 CQS 2 , 



-,1/2 



I 1/2 



xe -to' 1 

J 9 = p^ fed [co S (e/2) + B in(e/2)] 

2sin6» 



(8) 



For short wavelengths, kd ^> 1, the depen- 
dence r(k) is dominated by the integral I\ as 
-(X 2 /2)(kdshi8/ir)~ 3 / 2 e- 2kd / sind . In the long- 
wavelength limit, kd <C 1, one can replace the exponent 



-kd/ sin 9 



by 1. In what follows, we will focus on small 9, 



where d^ and r diverge. Note that the pole contribution 
in Eq. ([5]) diverges for 9 — > much stronger that the 
principal value contribution, which is oc ln(l/0), so 
that r !» ir 2 X 2 /29. We also notice that in the smail-6> 
domain, the relation r ~ dLdn holds. This relation can 
be understood from the following reasoning. 

There are two contributions to the reflected wave in the 
second Born approximation, (i) The wave deflected into 
the second tube to the right with the amplitude (solid ar- 
row in Fig. 1), undergoes a secondary deflection back into 
the first tube (dashed arrow in Fig. 1) with amplitude d^. 
(ii) The wave deflected into the second tube to the left, 
c?l , is subsequently deflected back into the first tube with 
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the amplitude da, Fig. 1. The sum of the two contribu- 
tions amounts to r = (ci + C2)di,dii. Remarkably, both 
numerical factors C\ and C2 are equal to 1/2. This is a 
consequence of a strong difference in distances at which 
formation of the primary left- and right-deflected waves 
takes place. The wave dh is formed within ~ 1/k from 
the intersection, while the wave, dfj, is formed within a 
much broader interval ~ l/(k8). Therefore, in second 
tube, at some distance y from the intersection, such that 
1/k <C y <C l/(k9), the amplitude of the left-deflected 
wave is already c£l, while the amplitude of the right- 
deflected wave is only \dn. Subsequent formation of the 
contribution {it) occurs at y ~ l/{k8), so that the corre- 
sponding amplitude is (|<ifi) dhi i- e -, c 2 = 1/2- On the 
other hand, formation of the contribution {i) takes place 
only over negative —l/{k6) < y < 0, and thus results in 
dh (\d R ), i.e., ci = 1/2. 

Semiclassical description. From Eq. ^ one can see that 
for 9 < 7rA the Born approximation renders an unphysi- 
cal result, namely, du > 1, suggesting that this approxi- 
mation is not applicable for small 9. This manifests the 
change in the mechanism of the plasmon scattering which 
takes place for 9 < A. Indeed, at small 9, incident wave 
travels closely to the wave da over a long distance, so 
that their amplitudes get mutually redistributed. Im- 
portantly, in describing this redistribution one can: (i) 
neglect both left-deflected, dr,, and reflected, r, waves 
and (ii) employ semiclassical approach, which yields 



'2A f°° 
t{k) = cos ( — / dzKo 



\Jk 2 d 2 + z 2 



(9) 



2A f°° r 

dn{k) = isin ( — J dzK Q \f k 2 d 2 



where Kq is the MacDonald function. A remarkable fea- 
ture of this result is that t and r oscillate strongly with 9, 
and that the oscillations scale with the interaction con- 
stant. Note, that in the short-wavelength limit kd 1, 
the perturbative result Eq. §§§ is valid even for A > 9. 
Using the large- argument asymptote of K${z), it is easy 
to see that Eq. © reproduces Eq. © in this limit. For 
long- wavelengths, kd <C 1, Eq. ([9]) yields d^ = sin(7rA/0), 
t = cos{tt\/9), so that the perturbative and semi-classical 
results match at A/0 < 1. 

To outline the derivation of Eq. © , we note that the 
system of equations, Eqs. ([2]), can be rewritten as two 
independent closed equations, 

(uj 2 /v 2 f ) u ± {x) = D{u ± {x)} ± F{u±{x)}, (10) 

where combinations u±(x) = u%{x) ± 1*2(2;) are intro- 
duced. Searching for the semiclassical solution of Eq. (fTU]) 
in the form u±{x) — e lkx + l v± ( kx ) f with slowly varying 
phase, ip' ± <C 1, we find 



2<p' ± {kx) =tA^ (l + <p' ± (kx)) ky/d 2 + x 2 9 2 . (11) 



In evaluating the r.h.s. we assumed that 9 is small. We 
see that when A is small, the assumption, <£>' ± -C 1, is 
justified. Then the smallness of ip' ± allows one to neglect 
it in the argument of K . Upon integrating Eq. (JTTJ, 
we find tp±. Then transforming back to u\ and 1*2, we 
recover Eq. The expression for r{k) generalized to 
the domain 9 < A < 1 follows from Eqs. ([6]) and (J9]) 



r{k)\ 



dnd L 



(12) 



_ 2kd 
e ~ B ~ sin 



— / dzKo 




\Jk 2 d 2 + z 2 



and in the long-wavelength limit simplifies to r\e < \ — 
(7rA/2)sin(7rA/6»). 

Tunnel density of states. Most importantly, the non- 
trivial dependence of the plasmon scattering on A and 
9 manifests itself in the observables, e.g., in the depen- 
dence of tunnel density of states, v(co, x), on the distance, 
x, from the intersection. To illustrate this, consider first 
a single nanotube with inhomogeneity at x = which 
scatters plasmons with reflection coefficient f{k). Then 
the correction to the tunnel density of states reads 



Sv{uj, x) 
v Q {u) 



= T{a + l)^Ja 2 +a/2 



(13) 



\r{u/s)\ 
{2uox/s) a+ 



2ujx 
s 



ip{w/s) 



7TQ 
~~2 



where f{k) = arg(f). Eq. (|13p follows from the expres- 
sion for interaction contribution to the local Green func- 
tion which takes into account the plasmon scattering, 



5{x,t) = 
exp < — 7r / dk 



(14) 



\J x u k \ 2 sk \u k {x)\ 2 v B 



8v f 



Here u k {x) are the plasmon eigenmodes: u k {x) = (e tkx + 
r{k)e- ikx )/V2n, and u k {x) = ^{1 - \f \ 2 )/2tt e lkx , for 
kx < and kx > 0, respectively. Expanding the ex- 
ponent in Eq. (|14|) with respect to f, and evaluating 
v{uj,x) = TT-^e J °° dte iuJt g{x,t), we arrive at Eq. (p]l. 
Eq. (fT4|) emerges upon representing electrons via the 



dual bosonic fields 9i a and 4>%a 1 ipi 



labels the two bands, a =f,l are the spins [8]. In- 
teraction is completely described by the charged field 
9 C = i J2ia ® ia ' wri ile the three neutral sectors are non- 
interacting. Expansion [l5| 9 c (x) = ir^/no J dku k {x)Q k , 
<t> c (x) = j-j^ J dk[J x dyu k (y)]P k over the plasmon eigen- 
modes, u k (x), reduces the interacting Hamiltonian to a 
system of harmonic oscillators {Q k , Pk} yielding Eq. (TIT 



A simple reasoning allows to generalize Eq. (TT3)) to the 
case of two intersecting nanotubes. Indeed, with inter- 
section playing the role of inhomogeneity, instead of one 
reflected wave with reflection coefficient r we have two 
independent modes, u±(x), solutions of Eq. (jlOp . with 
reflection coefficients r±. It is important that while the 
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absolute values r + and r_ are the same and equal to |g?l|, 
given by Eq. ([6j> , their phases are different and are equal 
to ip + (kx) — ip + (—kx) and (fi-(kx) — ip_(—kx) + n, re- 
spectively, where (f± are determined by Eq. (fTTj) . Due to 
this difference in phases, the oscillations oc sm[2ujx/ s + ip] 
in Eq. (fT"3")) transform into a beating pattern 



5v(lo, x) 
I 

x sin 



\ -2ud/sB 

= -T(a+l)^/^ + a/2— ^ ■ (15) 



(2lux/s) 



ujx9 I s 

— f dzKo 



V 



/ dzK [v/M/.s) 2 + ^ cos - f ) 



Eq. (|15|) is our main result. Remarkably, the shape of the 
envelope of cos(2wa;/s) oscillations depends strongly on 
the interaction parameter, A, offering a unique signature 
of Luttinger liquid behavior. In particular, the number 
of nodes in the envelope is equal to X/8. Examples of 
oscillations Eq. (fT"5)) are plotted in Fig. 2 in terms of tun- 
nelling conductance, G(V,x) oc v(w = V,x), for different 
interaction parameters. Note, that the language of re- 
flected plasmons, r+, r_, applies at distances x 3> s/V, 
over which the reflection coefficient is formed. Since the 
characteristic scale of the envelope is s/V 6, Eq. (| 1 5|) is 
valid as long as 9 <C 1. For large x 3> s/V '8, the ar- 
gument of the sine in Eq. (| 1 5|) saturates at ttX/9. On 
the physical grounds, the magnitude of the cos(2Vx/s) 
oscillations at large x should be given by Eq. (fl~3|) . with 
element of scattering matrix r instead of r. From Eq. (fl~2|) 
we realize that this is indeed the case. 
Implications. Our main prediction is that, for purely 
capacitive coupling between nanotubes, the conductance 
GiV) into one or both ends of crossed nanotubes must 
exhibit a structure, like shown in Fig. 2, with a large 
characteristic "period" V ~ s/(L6), where L is the dis- 
tance from the end to the intersection. Smallness of 9 
insures that this structure (envelope in Fig. 2) is distin- 
guishable from size-quantization-like "filling" of the en- 
velope [1, [3], which changes with the period V = irs/L. 
Also, an important prediction is that the envelope beat- 
ing structure in Fig. 2 vanishes with temperature much 
slower than the filling, which vanishes at T ~ s/L. 

The loop geometry of Ref. 17 offers another possible 
experimental implication. For this geometry, the easi- 
est way to compare the Sagnac oscillations in Ref. Qj] 
and our finding Eq. (|15p is to assume that interaction 
is weak. Then the contribution to the differential con- 
ductance from the Sagnac effect is roughly the product 
of "size-quantization" oscillations, oc cos(2VL/vp), and 
the envelope oc cos(2V Lu g /v F ), where V is the source- 
drain bias; L is the loop perimeter and u g is the gate- 
induced detuning of the "left" and "right" velocities. Our 
Eq. (fT5|) for this geometry contains the same first cosine 
cos(2V L/vf), while the envelope is entirely due to inter- 
actions. Thus, common feature of the two effects is that 



envelopes survive at "high" temperatures when Fabry- 
Perot oscillations vanish. 

Concluding remarks. Adding a second parallel nanotube 
to a given one leads [HI to a reduction of a in G(V) by 
a factor of 2. One could expect that for a finite crossing 
angle, the effect of the second nanotube is weaker. We 
found, however, that GiV) depends on 9 in a nonanalyt- 
ical fashion when 9^0. This nonanalyticity translates 
into a peculiar bias dependence of G, as shown in Fig. 
2. Thus, for crossed nanotubes, G(V) is extremely sensi- 
tive to the value of intratube Luttinger liquid parameter, 
g. In armchair nanotubes, the currently accepted value 
HSU is in the range 0.19 -j- 0.26. We emphasize that 
changing g from 0.19 to 0.26 leads to the increase of the 
interaction parameter, A, by a factor of 2, which would 
have a drastic effect on the shape of envelope in 8G(V), 
Fig. 2. 

Concern ing relation between our study and earlier 
studies [H, [HI of crossed nanotube junctions, this rela- 
tion is exactly the relation between plasmon and elec- 
tron scattering. In the above papers anomalies were due 
to either direct passage of electrons through the crossing 
point [l3T | or due to crossing- induced electron backscatter- 
ing ||. Scattering of plasmons was disregarded in Ref. 5. 
This is justified for perpendicular nanotubes of Ref. ]5. 
As shown in our manuscript, scattering of plasmons be- 
comes important at small angles. 
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